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(54) Method of compensation for the thickness of an organ 



(57) A method of compensation for the thickness of 
an organ in an X-ray machine of the type comprising an 
X-ray source and a means of detection of the beam of 
X-rays after it has passed through the said organ, the 
detection means being capable of converting the X-ray 
beam into a digital electronic signal, in which, from a 
digitized image, an image of the radiologic thicknesses 



of an organ traversed by the X-ray beam is calculated, 
a thickness threshold is defined, from this there is de- 
rived an algebraic compensation image to bring the pix- 
els of a level below or above the said threshold back to 
the value of the threshold, and the thickness image and 
a proportion of the compensation image are summed 
up, to obtain a compensated thickness image. 
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Description 

[0001] The present invention relates to the field of radiologic imagery which makes it possible to visualize an organ 
or an organ part, in general of the human body. 

s [0002] Conventionally, X-ray radiography has been performed by means of sensitized films acted upon by X-rays 
after their passage through the organ to be studied. Radiologists have been trained in the interpretation of such images. 
Novel imagery technologies - solid state detector and digital acquisition system - must be adapted to current practice 
and must furnish an equivalent perception of pertinent information that radiologists have to examine. In particular, one 
of the requirements that digital systems must satisfy consists in an extension of the dynamics of gray level that simulates 

10 a conventional film as faithfully as possible. To that end, the digital image is displayed on a screen that the radiologist 
adjusts in interactive fashion to identify all the clinical signs in perceiving relations among the various components of 
the image. But the precise perception of density information by means of the image displayed on the screen is limited 
by the dynamics of the screen. Now one must automatically pass from the high contrast of the image of origin, for 
example in the region of 30 to 50, to the low contrast' offered by a video screen. 

75 [0003] Currently, physicians apply conformation techniques to patients with the aim of limiting the dynamics of ac- 
quired images, for example by compression of thick regions, and/or by addition of absorbent substances to compensate 
for zones of little thickness. In the case of X-ray mammography, the breast is compressed to the smallest and most 
constant thickness possible. In the field of cardiology, contour filters are used to avoid problems related to poorly 
absorbent zones of the chest, such as the lungs. These filters consist in plates of a shape complementary to that of 

20 the heart, made of materials of given coefficients of X-ray absorption. 

[0004] However, these techniques prove to be inadequate for the acquisition and digital processing of images, and 
difficult to carry out. 

[0005] The present invention seeks to remedy the disadvantages described above. 
[0006] The present invention is to remedy the disadvantages described above. 
25 [0007] The present invention compensates for variations in thickness, in particular at the limit between regions of 
great density and regions of low density of the radiographed organ. 

[0008] . The method of compensation for the thickness of an organ is provided by an X-ray machine of the type con- 
taining a source of X-rays and a means of detection of the beam of X-rays after it has passed through the said organ, 
the means of detection being capable of converting the X-ray beam into a digital electronic signal. From a digitized 

30 image, there is calculated an image of the radiologic thicknesses of the organ traversed by the X-ray beam, the thickness 
zero corresponding to the zones of the image without organ, a thickness threshold is defined, from this there is derived 
an algebraic compensation image making it possible to bring the pixels of lower or higher level back to the threshold 
to the value of the threshold, and the thickness image and a proportion of the compensation image are summed up, 
to obtain a compensated thickness image. 

35 [0009] Advantageously, the compensation image is filtered through a low-pass filter, then through the mask, prior to 
summation of the images. 

[0010] Thus, the invention offers an image processing method furnishing an image equivalent to that which would 
be obtained by arranging an absorbent liquid on the edges of the organ over a part of their height and makes it possible 
to simulate a physical phenomenon which furnishes a better understanding of the optimal adjustment of the parameters 
40 that must be effected. This image processing gives the organ studied a natural and pleasant appearance. The method 
is easy to carry out by the adaptation of existing digital X-ray machines and may be applied to the radiography of any 
organ whatsoever. 

[0011] The invention will now be described in greater detail, by way of example, with reference to the drawings, in 
which:- 



45 



50 



55 



Figure 1 is a schematic representation of an organ composed of zones of various thicknesses; 

Figure 2 is a schematic representation of an organ composed of zones of various compensated thicknesses, 
according to the invention; 

Figure 3 is a schematic view of a breast in a compression plate; 

Figure 4 is a diagram representing the evolution of the thickness of the breast and the corresponding intensity of 
X-rays and corresponding to Figure 3; 

Figure 5, similar to Figure 3, shows the arrangement of an absorbent liquid in complementary shape with the breast; 
Figure 6 is a diagram similar to that of Figure 4 and corresponding to Figure 5; 
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Figure 7 is a diagram of the image processing; and 

Figure 8 is a histogram of thickness of the breast. 

5 [0012] The thickness of an organ as it is measured by X-rays is called radiologic thickness, in other words, taking 
account of the absorption of the materials traversed. For example, 1 cm of bone has the same radiologic thickness as 
4 cm of water. 

[0013] The thickness image may be obtained by means of Lambert's law: I = loe-u.t with I the number of photons 
received at a given point of the detection means, lo the number of photons not having passed through the organ and 
10 received at a given point of the detection means, u. the coefficient of linear absorption of the X-rays through the material 
traversed, and t the thickness of material traversed, from which there is derived: 1n lo - 1n I = ui; the product ui 
corresponding to the gray level for one pixel of a radiologic thickness image. 

[001 4] In one embodiment of the invention, a threshold is established corresponding to a given thickness depending 
upon the nature of the organ, all thicknesses above this threshold and then a predetermined part of the lesser thick- 
15 nesses remaining and closest to the threshold are suppressed, a mask defining the contour of the said organ is thus 
obtained, the value 1 is applied to the zones of the image corresponding to the said organ and the value 0 is applied 
to the zones of the image not corresponding to the said organ, and then the product of the mask by the digitized image 
is taken to obtain an image corresponding only to the organ. 

[0015] In one embodiment of the inventbn, a supplementary thickness threshold is defined as a function of a pre- 
20 determined proportion of pixels of a level below the said threshold, from this there is derived a compensation image 
making it possible to bring the pixels of a level below the said threshold back to the value of the said threshold, the 
compensation image is filtered through a low-pass filter and then through the mask, and the image corresponding only 
to the organ and a proportion of the compensation image are summed up. 

[0016] In one embodiment of the invention, the low-pass filter is calculated to eliminate the pixels corresponding to 
25 an organ filled with contrast medium in the compensation image and thus preserve them in the compensated image. 
[0017] Advantageously, erosion of the mask is effected prior to the product of the mask by the digitized image. 
[0018] In one embodiment of the invention, the predetermined part of the remaining thicknesses below and closest 
to the threshold is comprised between 1 and 5%. 

[0019] In one embodiment of the invention, the threshold is determined as a function of the size of the organ, for 

30 example 1 cm for an organ 5 cm thick. 

[0020] In one embodiment of the invention, the supplementary thickness threshold is defined as a function of a 
predetermined proportion of pixels equal to 20% of the pixels above the supplementary threshold. 
[0021] In another embodiment of the invention, the thickness of the organ is determined by comparison of the thick- 
ness histogram with a preestablished model curve, a parabola approaching the first of the maximums of the thickness 

35 histogram is calculated and the thickness threshold at the intersection of the parabola and the abscissa of the thick- 
nesses is selected. 

[0022] In one embodiment of the invention, two thickness thresholds, high and low, are defined, the compensation 
image being calculated to bring the pixels of a level below the threshold back to the value of the said low threshold, 
and the pixels of a level above the high threshold back to the value of the said high threshold. In the latter case, the 

40 compensation image has a negative value. 

[0023] As can be seen in Figure 1, the dynamics of an image of an organ may be greater than the dynamics of 
thicknesses representable by the visualization screen and which are represented by the two horizontal dotted lines. 
Thus, the blocks 1 to 3 are situated at too high a gray level while the blocks 8 and 9 are situated at too low a level. The 
physical phenomena or the parts of the organ that are located in such zones therefore are not visualized in satisfactory 

45 fashion. While preserving the same screen dynamics, first of all the blocks 5 to 9 could be visualized and then, sub- 
sequently and after modification of the parameters of luminosity of the screen, the blocks 1 to 4. Taking the example 
of mammography, a first image would be obtained on which only the contour of the breast would be seen, without its 
internal tissues being seen. The second image would show only its internal tissues, adipose zone or glandular zone, 
without its contour being apparent, which is hardly practical and complicates the radiologist's work. 

50 [0024] In the present invention, digital compensation of the raw image is based on the following principle: the method 
of image compensation simulates the addition or removal of an appropriate quantity of matter in border regions, for 
example between the zones of strong absorption and the zones of weak absorption of X-rays, so that the dynamics 
are reduced while taking account of the constraint of preservation of the differences and the real ratios among the 
anatomical structures. Thus, it can be seen that in Figure 2, blocks 1 to 9 are contained within the limits of the dynamics 

55 of the screen. Blocks 1 to 3, instead of being saturated in white, are brought back to just below the upper limit of 
saturation and preserve their original relation, namely: block 1 lighter than block 2 lighter than block 3, etc. The same 
is true for blocks 8 and 9, which, instead of being black, are situated near the lower limit of the dynamics, so as to be 
visible and while preserving their original difference, block 7 lighter than block 8 lighter than block 9. 
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[0025] Figure 3 shows a breast 1 compressed between a plate 2 and table 3 which make up part of an X-ray machine, 
not represented. The breast is subjected to an X-ray beam 4 emitted by a source not represented. A digital detector, 
not represented, is arranged on the path of the X-ray beam 4 after its passage through the breast 1 . It is noted that 
certain X-rays, labelled 5, pass through only a part of the thickness of the breast because of its rounded shape at the 
s end. The result of this is that these X-rays are less attenuated than those that have passed through the totality of the 
thickness of the breast, which risks leading to saturation of the image portion corresponding to the rays 5 and to display 
of a portion of black screen. In fact, a white portion corresponds to a thick zone and a black portion corresponds to a 
zone that js not very thick, by convention. 

[0026] The diagram of Figure 4 illustrates this situation. When the thickness diminishes, the intensity (or number of 
10 photons) increases considerably, which makes the edge of the breast difficult to see. The result of this is that tumors 
and microcalcifications at the edge of the breast may be missed by the radiologist analyzing- the anatomy of the central 
part of the breast. The method of compensation according to the invention makes it possible to bring the visibility of 
the edge of the breast to almost the same level as the visibility of the central part of the breast. 
[0027] To that end, fill-in is simulated to a level H (Figure 5) of a sector previously filled with air in contact with the 
is end of the breast. Fill-in is effected with a solution 6 that increases the apparent thickness of the breast at its edge. 
This is equivalent to positioning contour filters or placing sacks of flour along the arms or neck of the patient for arte- 
riography to avoid the saturation of image in these regions of little thickness. Providing a liquid 6 of course permits 
perfect adaptation to the shape of the region that must be compensated in thickness. An apparent thickness and an 
intensity whose evolution is represented in Figure 6 and which presents a much more favorable configuration are thus 
20 obtained. However, providing a liquid imposes handling constraints that it would be desirable to avoid. The method 
according to the invention makes it possible to simulate the presence of the liquid 6. 

[0028] The various steps of the method according to an embodiment of the invention are represented in Figure 7. 
First, the input image coming from the detection means is subsampled, for example by a factor of Z=4, to reduce 
execution time, see block 7. In one variant, the method may be carried out with the original resolution of the image to 

25 avoid any artifact in the appearance of the skin of the organ. 

[0029] In block 8, the logarithmic function that makes it possible to pass intensities is applied to the radiologic thick- 
nesses. In block 9, extraction of the breast is effected. In block 10, estimate of the minimum thickness of the adipose 
zone of the breast is effected. In block 11 , the compensation function is calculated. In block 12, the image is reestab- 
lished in its resolution of origin. In block 1 3, the exponential function is applied. 

30 [0030] For a given point of the image, the number I of photons received by a given point of the detection means is 
determined by the following equation: I = I max x e-u.t, t being the thickness of the organ traversed by the X-rays and 
u. being the density of the organ traversed. Application of the logarithmic function results in: u. x t = 1n (Imax) - 1n (I). 
Imax is known and corresponds to the number of photons that are received on a given point of the detection means 
and have not passed through the organ. Thus, the product u. x t, which is called densitometric or radiologic thickness, 

35 is found. 

[0031] In practice, it may be preferable to apply a slightly different equation based on the magnitude G, which is the 
gray level of a given point of the detection means: u.t = K x 1n (Gmax + 1) - K x 1n (G + 1). The gain K is applied to 
protect the variations in gray level greater than 1% of the maximum intensity, such as K x 1n (Gmax + 1) - K x 1n (0.99 
Gmax + 1) £ 1. Application of this logarithmic function makes it possible to obtain a histogram of the thickness, as 
40 illustrated in Figure 8. 

[0032] Then the minimum thickness of the adipose zone is estimated. The histogram of thickness of the breast 
(Figure 8) includes three principal regions, the base, the skin or edge of the breast and the body of the breast proper. 
First of all, a thickness threshold is defined, preferably adjustable as a function of the dimensions of the organ, for 
example equal to 1 cm for mammography. All thicknesses greater than this value are eliminated, then a given percent- 

45 age, for example 1 %, of thicknesses below the threshold value and closest to the said threshold value are eliminated. 
Then, a coefficient 1 is assigned to the zones corresponding to the breast, i.e., the zones previously excluded, and a 
coefficient 0 to the zones not corresponding to the breast, i.e., the zones previously preserved. A mask is thus obtained 
which, by multiplication with the image of the breast, makes it possible to preserve only the pixels corresponding to 
the breast and to suppress ail the others. 

50 [0033] Then, an erosion (operation of mathematical morphology) is carried out on a given thickness, for example 
about 1 cm thickness in the case of mammography, of this mask. The new mask thus obtained corresponds to the 
body of the breast. Then the product of this mask by the input image is taken, which makes it possible to preserve the 
image of only the region of the body of the breast, then 20% of the weakest pixels are eliminated, which corresponds, 
in Figure 8, to suppressing 20% of the surface of the thickness histogram in the region of the body of the breast closest 

55 to the skin region. The limit of this 20% corresponds to the value of minimum thickness of the adipose zone ua x ta min. 
[0034] Then the compensation to be provided is calculated. First of all, the image previously calculated is passed 
m through a low-pass filter to preserve only slow variations of thickness u, x t! The compensation image uaH is defined 
as being that which must be added to the filtered image u. x t to obtain values equal to the minimum thickness of the 
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adipose zone uata min: uaxH = uaxtamin-u.xt>Q. The liquid 6 whose presence is simulated to a height H has 
the same density na as the adipose zone. 

[0035] Multiplication of the compensation function by the mask is effected before and after low-pass filtering in order 
to avoid outside-the-breast compensation. The image u. x t and a proportion of the compensation image uaH are 
s summed up in order to preserve the hierarchy of the compensated thicknesses. In other words, compensation is limited 
to a rate of less than 100%. 

[0036] Of course, the parameters described above may be subjected to various adjustments in order to adapt them 
to various adipose zones of the breast or to other organs. 

[0037] In a variant, use of the mask may be replaced by a simulation using an organ model. Then the thickness of 
10 the organ is determined by comparison of the thickness histogram with the preestablished model curve, a parabola 
approaching a part of the thickness histogram is calculated and the thickness threshold at the intersection of the pa- 
rabola and the abscissa of thicknesses is selected. 

[0038] In a variant, taking the product of the mask by the compensated image instead of applying it to the compen- 
sation image could be considered. This would have the drawback of suppressing certain variations of the base, as well 

is as effecting compensation of parts of no interest. 

[0039] In the field of vascular imagery, in particular for examinations of cardiac angiography, the same method may 
be applied by compensating poorly absorbent zones of the chest by virtual addition of matter and by virtual removal 
of matter in the dense zone of the diaphragm. In similar fashion, operations of arteriography of the carotid or of the 
lower extremities may advantageously be performed with the present invention. Then a high threshold and a lowthresh- 

20 old are defined, the compensation image being calculated to bring the pixels of a level below the low threshold back 
to the value of the said low threshold, and the pixels of a level above the high threshold back to the value of the said 
high threshold. 

[0040] For thick zones such as the chest, the following equations may be used: 

25 

G = Gmax x e-u.t 



|it = 1n(Gmax)-1n(G) 

30 

u.aH = compensation in the lung 



^ fibL = compensation of the diaphragm or of the vertebrae 

[0041] Compensation is effected according to the following rule by adding the compensation uaH and subtracting 
the compensation ubL: 

40 

u.t + jxaH jibL, 

which results, after application of the exponential function, in 

45 

G x e-uaH x e+u.bl_ 

[0042] Thus, in the dense zone of the diaphragm, compensation carried out according to the invention will correspond 
to a virtual removal of matter so as to bring the dynamics of the image back to a value corresponding to that of the 

so visualization screen. The same will be true in the zone of low density of the lung, where matter will be virtually added. 
[0043] To preserve satisfactory visualization of arteries that are injected with contrast liquid but are relatively small 
in diameter, the filter will be adjusted so that it does not see the arteries, which wilt thus be absent from the compensation 
image and preserved in the compensated image. For example, if the filter eliminates details smaller than or equal to 
1 cm, these will be preserved in the compensated image, while the details of dimensions greater than 1 cm will be 

S5 erased. 

[0044] The invention makes available a method of image processing capable of simulating a physical operation of 
provision or removal of an absorbent liquid, which permits ready comprehension of its functioning. The image obtained 
according to the method exhibits no artifacts and is natural in appearance. The method may be carried out in various 
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fields of radiology, while eliminating the need for physical compensation devices. 



Claims 

5 

1. A method of compensation for the thickness of an organ in an Xray machine, of the type comprising a source of 
X-rays and a means of detection of the X-ray beam after it has passed through the said organ, the detection means 
being capable of converting the X-ray beam into a digital electronic signal in which, from a digitized image there 
is calculated an image of the radiologic thicknesses of the organ traversed by the X-ray beam, comprising the 
10 steps of: 

defining thickness threshold; 

deriving an algebraic compensation image to bring the pixels of lower or higher level back to the said threshold 
is to the value of the threshold; 

summing the thickness image and a proportion of the compensation image to obtain a compensated thickness 
image. 

20 2. The method according to claim 1 , comprising the steps of: 

establishing a threshold corresponding to a given thickness depending upon the nature of the organ; 

suppressing all thickness above this threshold and then a predetermined part of the lesser thicknesses re- 
25 maining and closest to the threshold to obtain a mask defining the contour of the said organ; 

applying, the value 1 to the zones of the image corresponding to the organ and the value 0 to the zones of 
the image not corresponding to the organ; and 

30 taking the product of the mask by the digitized image to obtain an image corresponding only to the organ. 

3. The method according to claim 2, in which a supplementary thickness threshold is defined as a function of a 
predetermined proportion of pixels of a level above the threshold, from this there is derived a compensation image 
to bring the pixels of a level above the threshold down to the value of the threshold, the compensation image is 

35 filtered through a low-pass filter and then through the mask, and the image corresponding only to the organ and 

a proportion of the compensation image are summed up. 

4. The method according to claim 3, in which the low-pass filter is calculated to suppress the pixels corresponding 
to an organ filled with contrast medium. 

40 

5. The method according to any one of claims 2 to 4, in which an erosion of the mask is performed prior to the product 
of the mask by the digitized image. 

6. The method according to any one of claims 2 to 5, in which the predetermined part of the remaining thicknesses 
45 below and closest to the threshold is equal to 99%. 

7. The method according to any one of claims 2 to 6, in which the threshold is determined as a function of the size 
of the organ. 

so 8. The method according to any one of the preceding claims, in which the thickness threshold is defined as a function 
of a predetermined proportion of pixels equal to 20%. 

9. The method according to claim 1 , in which the thickness of the organ is determined by comparison of the thickness 
histogram with a preestablished model curve, a parabola approaching the first maximum of the thickness histogram 

55 js calculated and the thickness threshold at the intersection of the parabola and the abscissa of the thicknesses 

is selected. 

10. The method according to claim 1 , in which two thickness thresholds, high and low, are defined, the compensation 
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image being calculated to bring the pixels of a level below the low threshold back to the value of the said low 
threshold, and the pixels of a level above the high threshold back to the value of the said high threshold. 

11. The method according to any one of the preceding claims, in which the compensation image is filtered through a 
low-pass filter, then through the mask, prior to summation of the image corresponding only to the organ and the 
compensation image. 
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FIG.7 




11 



EP 1 004 985 A1 



FIG.8 
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